Further study of double charmonium production in e^e annihilation at Belle 
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We report a new analysis of double charmonium production in e^e^ annihilation using a data 
sample collected by the Belle experiment. We confirm our previous observation of the processes 
e"'"e~ J/ipric{xcQ,Vc{'^S)) and perform an angular analysis for these processes. Processes of the 
type e'^e~ —^ tj){2S){cc)res are observed for the first time. We also observe a new charmonium state 
- X(3940), produced in the process e+e" J/i/;X(3940). 
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PACS numbers: 13.66.Bc,12.38.Bx,14.40.Gx 



The surprisingly large rate for processes of the type 
e+e^ ^ J/tjjrjc and J/^ (cc)non-rcs observed by Belle [ij 
remains unexplained. In the Belle analysis with a data 
sample of 45 fb"-'^, the presence of the process e+e^ 
J l^p rjc was inferred from the rj^ peak in the mass spec- 
trum of the system recoiling against the reconstructed 
J/^ in inclusive e^e^ J/ipX events. Following the 
publication of this result, the cross-section for e+e^ — > 
J/'0 Tjc via e^e~ annihilation into a single virtual photon 
was calculated to be ~ 2 fb Q , which is at least an order 
of magnitude smaller than the measured value. Several 
hypotheses have been suggested in order to resolve this 
discrepancy. In particular, the authors of Ref. |^ have 
proposed that processes proceeding via two virtual pho- 
tons may be important. Other authors 01 suggest that 
the final states observed by Belle contain a charmonium 
state and a M ~ 3GeV/c^ glueball, which are predicted 
by lattice QCD. Possible glueball contributions to the 
Xco signal are also discussed in Ref. 5] . In this paper we 
report an extended analysis of the e'^e" — > J/'0(cc)rcs 
process to check the above hypotheses and provide extra 
information that might be useful to resolve the puzzle. 
This study is performed using a data sample of 155 fb~^ 
collected around the T(4S') resonance with the Belle de- 



tector at the KEKB asymmetric energy 
rings. 
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The analysis procedure is 
Refs. [1 The recoil mass M, 
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J/tp center-of-mass (CM) energy and momentum, respec- 
tively. The MiccoiK^/V') spectrum for the data is pre- 
sented in Fig. ^ clear peaks around the nominal rjc and 
XcO masses are evident; another significant peak around 
~ 3.63GeV/c^ is identified as the j]c{'2S). The authors 
of Ref. estimated that the two-photon-mediated pro- 
cess e^e~ J/ipJ/il) has a significant cross-section. To 
allow for a possible contribution from the exchange of 
two virtual photons, we fit the spectrum in Fig.^includ- 
ing all of the known narrow charmonium states below 
DD threshold. The fit results are listed in Tabled The 
yields for r^c, XcO, and r]c{2S) have statistical significances 




i(J/v) 



3.8 

(GeV/c^) 



FIG. 1: The mass of the system recoiling against the recon- 
structed J/ip. The curves are described in the text. 



TABLE I: Summary of the signal yields (N), charmonium 
masses (A/) and significances (it). 



(cc)ros 


N M[GeV/c^] 


a 


Vc 


235 ± 26 2.972 ± 0.007 10.7 


J/i> 


-14 ±20 fixed 




XcQ 


89 ± 24 3.407 ± 0.011 


3.8 


Xcl+Xc2 


10 ± 27 fixed 






164 ± 30 3.630 ± 0.008 


6.0 


^{2S) 


-26 ± 29 fixed 





between 3.8 and 10.7. The fit returns negative yields for 
the J/'0 and ■0(25'); the Xci and Xc2 yields are found to 
be consistent with zero. A fit with all these contributions 
fixed at zero is shown as a solid line in Fig.^ the dashed 
line in the figure corresponds to the case where the con- 
tributions of the J/ij}, Xci, Xc2 and ^{2S) are set at their 
90% confidence level (CL) upper limit values; the dotted 
line is the background function. 

Given the arguments in Ref. it is important to 
check for any momentum scale bias that may shift the 
recoil mass values and confuse the interpretation of peaks 
in the Mj-ocoii spectrum. We use e+e" V(2S')7, 
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TABLE II: The a parameters obtained from fits to the 
production and helicity angle distributions for e+e^ 

J/lp (cc)res. 



(cc)rcs 


Separate 

Qprod Ohcl 


Simultan. 

aiicl = Qprod 


XcO 

??c(2S') 


1 4+11 Q t-+U.7 
i-4_0.8 U-O-O.S 
1 7+0.5 n 7+0.7 
i-'_0.5 U.(_Qg 

^■^-1.2 "■'^-0.7 


0.87l°;f3 



tp{2S) — > J/'07r+7r~ events to calibrate and verify the 
recoil mass scale. From the study of the spectrum of 
recoil masses squared against ipC^S) in the data, we cal- 
culate that the J/^p recoil mass is shifted by not more 
than 3MeV/c^. As an additional cross-check we fully 
reconstruct double charmonium events. We find 3 pure 
events with J/ip rjc combinations with energies consistent 
with total CM energy. Based on the rj^ yield in the 
Mrccoii{J / 4') distribution, we expect 2.6 ± 0.8 fully re- 
constructed events, consistent with the observed signal. 
Thus we conclude that the peak in Mrecoii(>//V') is domi- 
nated by J/ip rjc production. We also search for fully re- 
constructed events with J/^ J combinations and find 
no such candidates in our data. Based on the calibra- 
tion of the Mrccoi\{J / ip) scale, the result of the fit to 
the Mj-ccoii{J / "tp) distribution and the full reconstruction 
cross-check, we confirm our published observation of the 
process e+e^ J/tprjc and rule out the suggestion of 
Ref. 'i^ that a significant fraction of the inferred J/^ rjc 
signal might be due to J /ip J jip events. 

The reconstruction efficiencies for the J /ipric, J/'ipXcO, 
and J/tpr]{2S) final states strongly depend on Op^od, the 
production angle of the J/tp in the CM frame with re- 
spect to the beam axis, and the helicity angle ^hei- We 
therefore perform an angular analysis for these modes 
before computing cross-sections. We fit the MrGcoii( J/V") 
distributions in bins of |cos(6'prod)| and |cos(0hci)|, and 
correct the yield for the reconstruction efficiencies deter- 
mined bin- by-bin from the MC. The results are plotted 
in Fig. 13 together with fits to functions A{1 + a cos^ 9) 
(solid lines). We also perform simultaneous fits to the 
production and helicity angle distributions for each of 
the (cc)res states, assuming J/ip (cc)rcs production via a 
single virtual photon and angular momentum conserva- 
tion, thus setting aprod = ahci- The values of the param- 
eter a are listed in Tabled The angular distributions 
for the J/iprjc and J/Tpr]c{2S) peaks are consistent with 
the expectations for production of these final states via a 
single virtual photon, aprod ~ cthci ~ +1 The predic- 
tion for a spin-0 glueball contribution (e+e~ — > J/ipQo) 
to the J/t/jVc peak, aprod = cthoi — —0.87 0, is disfa- 
vored. The process e+e^ — » 7* — ^ J/ipXcO can proceed 
via both S- and D-wave amplitudes, and predictions for 
the resulting angular distributions are therefore model 
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FIG. 2: Distributions of cosines of the production (left) and 
J/ip helicity angles (right) for e+e~ — > J/tpric (top row), 
e^e" — > J/ipXcO (middle row) e'*'e~ — > J/ipr]{2S) (bottom 
row). The solid lines are results of the individual fits; the 
dotted lines are the simultaneous fit results. 

dependent. Our results disfavor the NRQCD expecta- 
tion aprod = tthei — 0.25 H, |^ , and are more consistent 
with 5- wave production, where aprod = OLhc\ = ^ 1- 

To calculate the cross-sections we fix the production 
and helicity angle distributions in the MC to 1 -I- cos^ 6 
for J / ip r]c{rjc{2S)) , and to 1 — cos^ 9 for J/ipXcO- To 
reduce the model dependence of our results due to the 
effect of initial state radiation, whose form-factor depen- 
dence on of the virtual photon is unknown, we cal- 
culate cross-sections in the Born approximation. As in 
Ref. 0, we present our result in terms of the product 
of the cross-section and the branching fraction of the re- 
coil charmonium state into more than 2 charged tracks: 
(7 X S>2, where S>2((cc)res) = S((cc)rcs ~* > 2 charged). 
The cross-sections are given in Table HVl 

We perform a similar study with reconstructed 
ip{2S) J/ipTT~^Tr~ decays to search for e+e^ — > 
ip{2S){cc)rcs processes. The recoil mass spectrum for the 
data is presented in Fig.|3 peaks corresponding to the rjc, 
XcO, and r]c{2S) can be seen. The fit to the Mrecoii(V'(2S')) 
distribution is identical to the Mrecoii("//V') fit, but due 
to the limited sample in this case, the masses of the 
established charmonium states are fixed to their nom- 
inal values; the ric{2S) mass is fixed to 3.630 GeV/c^ 
as found from the Afrccoii(<//'0) fit. The signal yields 
are listed in Table IIIII Significances for the individ- 
ual 77c, XcO, and rjc{2S) peaks are in the range 3 ~ 4(t; 
the significance for e+e~ ip(2S) (cc)res, where (cc)res 
is a sum over rjc, XcO, and ric(2S), is estimated to be 
5.3(7. In Fig. 121 the result of a fit with only 77c, XcO and 
ric{2S) contributions included is shown as a solid line; 
the dashed line shows the case where the J/i^, Xci, Xc2, 
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FIG. 3: The mass of the system recoihng against the recon- 
structed tp{2S). The curves are described in the text. 



TABLE III: Summary of the signal yields (A'') and signifi- 
cances (cr). 



(cc)rcs 


N 


a 




36.7 ± 10.4 


4.2 


J/V 


6.9 ±8.9 






35.4 ± 10.7 


3.5 




6.6 ±8.0 




?7c(25) 


36.0 ± 11.4 


3.4 


V'(25) 


-8.3 ±8.5 





and ip{2S) contributions are set at their 90% CL upper 
limit values; the dotted line is the background function. 
Finally, the calculated products of the Born cross-section 
and the branching fraction of the recoiling charmonium 
state into two or more charged tracks {a x S>o, where 
B>o((cc)ios) = B{{cc)^cs charged)) are presented 

in Table El 

Using even larger data set of 280 fb^^, which became 
available by the summer 2004, we extend the analysis 
of the recoil masses against J/tp above DD threshold. 
We find another significant peak around the mass of 



TABLE IV: Summary of the cross-sections for e"'"e~ ^ 
J/tp (cc)res and e+e~ — > ip{2S) (cc)rcs. B>2{>o) denotes the 
branching fraction for final states containing more than 2 (at 
least one) charged tracks. The units are fb, and the upper 
limits are set at 90% CL. 

(cc)rcs O-Born X S>2 UBorn X B>0 

~rj'c 25.6 ± 2.8 ± 3.4 16.3 ± 4.6 ± 3.9 

J/i^ < 9.1 < 16.9 

Xco 6.4 ± 1.7 ± 1.0 12.5 ± 3.8 ± 3.1 

Xci+Xc2 < 5.3 < 8.6 

Vc{2S) 16.5 ±3.0 ±2.4 16.0 ± 5.1 ± 3.8 

ip{2S) < 13.3 < 5.2 




FIG. 4: The mass of the system recoiling against the recon- 
structed J/ip- The curves are described in the text. 



M - 3.940 GeV/c^ (Fig. E)). We denote this peak as 
X(3940). The fit to this spectrum includes the signal 
function, which is a sum over four charmonium states: 
Vc: XcO: Vc{2S), X(3940), and the background function 
that includes also the possible contribution from e+e^ — > 
J/tpDD events. The mass of the new charmonium state 
found by the fit is Mx = (3.940 ±0.012) GeV/c^; the sig- 
nal yield is = 149 ± 33 events and the significance of 
the signal is 4.5 a. The intrinsic width of the state is con- 
sistent with zero within a large error due to poor Mj-ocoii 
resolution. We set an upper limit on F to be 96MeV/c^ 
at 90% CL. 

We also search for the decay of AT (3940) into DD and 
DD* final states by reconstructing one of D mesons and 
requiring the second D or D* in the recoil mass spectrum 
against reconstructed J/tpD combinations. We find a 
significant signal of e+e^ J/'0X(394O) process when 
the event is tagged as AT (3940) — > DD*, while no signal 
is found for X(3940) -> D'D. We thus conclude that the 
dominant A" (3940) decay mode is DD* , and this state has 
a different nature from the recently found enhancement 
in J/ipw mass distribution around the same mass . 

In summary, using a data set of 155 fb~^ we confirm 
our published observation of e^e^ — > J/tprjc, J/ipXcO 
and J/tp 1^^(23) and find no evidence for the process 
e+e^ — * J/i/) J/i/i. We have calculated the cross-sections 
for e+e^ J/^rjc, J/i'XcO, and J/ijjr]c{2S) with better 
statistical accuracy and reduced systematic errors and 
set an upper limit for ^(e+e" ^ J/ip J/^)x B{J/ip > 
2 charged) of 9.1 fb at the 90% CL. Although this hmit 
is not inconsistent with the prediction for the J/^ J/ ip 
rate given in Ref. (3j , the suggestion that a large fraction 
of the inferred J/ ^ r]c signal consists of J/V' J/'0 events 
is ruled out. We have measured the production and he- 
licity angle distributions for e+e" — > J/iprjc, J/ipXcOi 
and J/'ipr]c{2S); the distributions are consistent with ex- 
pectations for these states, and disfavor a spin-0 glueball 
contribution to the rjc peak. We observe ^(2S')(cc)rcs 
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production for the first time, and find that the produc- 
tion rates for these final states are of the same magnitude 
as the corresponding rates for J/'0(cc)res- FinaUy, using 
a larger data set we observe the new narrow charmonium 
state at the mass Mx = (3.940 ± 0.012) GeV/c^. 
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